High Resolution mid-Infrared Imaging of SN 1987A^ 



Patrice Bouchet^'^, James M. De Buizer^'^, Nicholas B. Suntzeff^, I. John Danziger^, 
Thomas L. Haj^vard^, Charles M. Telesco^ and Christopher Packham^ 



ABSTRACT 

Using the Thermal-Region Camera and Spectrograph (T-ReCS) attached to 
the Gemini South 8m telescope, we have detected and resolved lO/xm emission 
at the position of the inner equatorial ring (ER) of supernova SN 1987A at 
day 6067. "Hot spots" similar to those found in the optical and near-IR are 
clearly present. The morphology of the lO/im emission is globally similar to the 
morphology at other wavelengths from X-rays to radio. The observed mid-IR 
flux in the region of SN1987A is probably dominated by emission from dust in 
the ER. We have also detected the ER at 20/im at a Aa level. Assuming that 
thermal dust radiation is the origin of the mid-IR emission, we derive a dust 
temperature of 180^'^^Qi^, and a dust mass of 1. — 8. x lO'^M© for the ER. Our 
observations also show a weak detection of the central ejecta at lO/xm. We show 
that previous bolometric flux estimates (through day 2100) were not significantly 
contaminated by this newly discovered emission from the ER. If we assume that 
the energy input comes from radioactive decays only, our measurements together 
with the current theoretical models set a temperature of 90 < T < 100 K and 
a mass range of 10^"^ — 2. x IO^'^Mq for the dust in the ejecta. With such dust 
temperatures the estimated thermal emission is 9(±3) x 10^^ erg s~^ from the 
inner ring, and 1.5(±0.5) x 10'^^ erg s~^ from the ejecta. Finally, using SN 1987A 
as a template, we discuss the possible role of supernovae as major sources of dust 
in the Universe. 
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1. Introduction 

The circumstellar envelope (CSE) surrounding SN 1987A consists of an inner equatorial 
ring (ER) flanked by two outer rings (Burrows et al. 1995). These rings, created ~20,000 
years before the explosion of the supernova, form the waist and caps of an hourglass or 
"bipolar" nebula ~1 pc across, enclosing an HII region and expanding into a diffuse medium 
that terminates in a dense shell 4 pc in radius (Crotts & Heathcote 2000). 

The collision between the ejecta of SN 1987A and the ER predicted to occur sometime 
in the interval 1995-2007 (Gaensler et al. 1997; Borkowski et al. 1997) is now underway. 
"Hot spots" have appeared inside the ER (Pun et al. 1997), and their brightness varies on 
time scales of a few months (Lawrence et al. 2000). In the next few years, new hot spots 
will continue to appear as the whole inner rim of the ER lights up, probably brightening by 
a factor 1000 during the next several years (Luo et al. 1994). 

The spectral energy distribution longward of 5 fxm has evolved continuously in time. 
After day 300 a cold dust-like component appeared (Suntzeff & Bouchet 1990; Bouchet & 
Danziger 1993; Wooden et al. 1993). An asymmetry in the profiles of optical emission lines 
that appeared at day 530 showed definitely that dust had condensed in the metal rich ejecta 
of the supernova (Lucy et al. 1991). Although it was discovered via spectroscopy (Danziger 
et al. 1989), the presence of the dust could be easily inferred from the spectral energy 
distribution: as the dust thermalized the energy output, after day 1000, SN 1987A radiated 
mainly in the mid infrared (Bouchet et al. 1991). Unfortunately, after approximately day 
2000, the flux emitted from SN 1987A in the thermal spectral region was too weak to be 
detected by the instruments and telescopes available. SN 1987A was observed at day 4100 
with ISOCAM onboard the ESA ISO satellite (Fischera & Tuffs 2000). Bouchet (2003) 
reported a weak detection of the supernova at day 4300 with OSCIR at the CTIO 4-m 
telescope. Except for these two observations, there have been no other detections of the 
mid-IR emission from the ejecta/ring region of SN 1987A over the last flve years. 

There exist very few mid-infrared observations of supernovae in general. Therefore SN 
1987A, the closest known supernova in 400 years, gives us an opportunity to explore the 
mid-IR properties of supernovae and their dust with the help of the newest generation of 
large-aperture telescopes and sensitive mid-IR instrumentation. 



^ Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc. under cooperative agreement with the NSF on behalf 
of the Gemini partnership: the National Science Foundation (United States), the Particle Physics and As- 
tronomy research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), 
the Australian Research Council (Australia), CNPq (Brazil), and CONICET (Argentina). 
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2. Observations 

The new T-ReCS mid-IR imager/spectrometer at the Gemini-S 8m telescope offers a 
combined telescope and instrument with diffraction limited imaging ( 0.3" resolution) and 
superbly low thermal emissivity. On 2003 Oct 4 (day 6067), we imaged SN 1987A with 
T-ReCS as part of the instrument's System Verification program. In a 23 minute on-source 
co-added image in the AT-band filter (A 7.70-12.97 //m), we were easily able to detect and 
resolve the ER (Fig. 1). This image shows several luminous "hot" spots distributed over the 
ring. The calibrated flux density integrated within an aperture of 1.3" radius is 9.9 ± 1.5 
mJy. This absolute calibration has been made using X^ff = 10.36 //m. No color correction 
was applied and this would most likely increase the flux density. The standard star used for 
this and all calibrations was HD 37160, whose flux density was taken to be 8.77 Jy in the 
A-band. On 2003 Dec 1 (day 6125) a short exposure of the supernova in the Qa-band (A 
17.57-19.08 yum) led to a 3.7a detection of the resolved ER with a flux density of 50.6 ± 6.6 
mJy (this image will be published when reinforced by improved S/N ratio images). We used 
a CMa with a flux density of 44.3 Jy at 18.30 //m for the flux cahbration of that observation. 

In Fig. 2 we compare our 10/im data with data obtained in the Hel line at CTIO 
(Bouchet 2003) at day 5749 (Fig. 2a), and in the F656N fllter as observed by HST at day 
5555 (Fig. 2b; ESO/ST-ECF Science Archive Facility). This fllter includes the Ha and [N 
II] (A 6583 A) lines. Our mid-IR image is also compared to the Chandra Observatory X-ray 
image obtained at day 5791 (Fig. 3a) (Park et al. 2004b), and the Australian Telescope 
National Facility (ATNF) 18.5 GHz (12 mm) maximum entropy restored image at day 6002 
(Manchester et al. 2003) (Fig. 3b). There is good agreement in shape and size between our 
mid-IR image and images obtained at other wavelengths. The mean radii and approximate 
surface brightness distribution (brighter on the east side) of the ring are similar from the X- 
ray to the radio, demonstrating that the dust is co-extensive with the gas components. The 
origin of that brightness asymmetry may be related either to the asymmetric distribution 
of the ejecta and/or to the density variation in the CSM (Park et al. 2004a,b), or to a 
time-dependent effect caused by the tilt of the ER as discussed by Panagia et al. (1991). 

The most hkely source of mid-IR radiation is thermal emission from cool dust (see 
discussion below), whereas the X-ray radiation is thermal emission from very hot gas (Park 
et al. 2003, 2004a). The radio emission is hkely to be synchrotron radiation as is the case 
for Gas A (Dunne et al. 2003), who claim that if the dust were at temperature 100 — 150i^ 

it could be heated by collisions with fast-moving electrons and ions in the hot gas seen with 
Chandra. Park et al. (2004a) argue that, until 2002 May 15, hard X-ray and radio emissions 
were produced by fast shocks in the GS HII region while the optical and soft X-ray emissions 
came from slower shocks in the denser ER. Park et al. (2003, 2004b) report morphological 
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changes in their last Chandra image, which is an indication of the blast wave now appearing 

to encounter the ER in the western side a few years after it reached the eastern side, and 
note that as of 2002 Dec. 31 (day 5791) correlations between the X-ray and the optical/radio 
images are more complex than the above simple picture. 



3. Discussion 

The origin of the 10 fj,m emission for both the ER and the ejecta, may be (1) hne 
emission from atomic species, (2) synchrotron or free-free emission, or (3) thermal emission 
from dust. As for the ER, its detection at 20/im seems to suggest that the emission is thermal 
emission from dust. A black body of temperature T^ust = 1801^q can be fitted to our data. 
An order of magnitude estimate can be obtained for the dust mass from the formula: 

- k{X) X n X B^{\,T) 

where D is the distance to the supernova, Fi,{X) the observed flux at lO^m, k{X) the dust 
mass absorption coefficient at this wavelength, and tt x B^{X,T) the Planck function. The 
approximate range for graphite and silicate grains for k at lOyum is 100 — 2000 cm^ g~^. We 
use D — blA Kpc (Panagia, N. 1999). Considering the uncertainty in the temperature we 
derive from this formula a dust mass range of 1 — 80 x lO^^M©. Park et al. (2004b) reported 
a plasma temperature of T ~ 2.8 x 10^ K and a gas density n — 235 cm~^ for the hard 
component of their two-shock model. According to Dwek (1987) our inferred temperature 
indicates a high gas density (n > 100 cm^"^) and a large grain size in the ER, thus in good 
agreement with the Chandra X-ray observations. 

Possible scenarios for the latc-timc mid-IR emission in Type 11 supcrnovac (Graham 
& Meikle 1986; Gerardy et al. 2000) are (1) dusty ejecta, (2) an infrared echo, or (3) dust 
heated from circumstellar interaction. In the last case, the dust could be (i) pre-existing 
dust in the GSM heated by the outer blast wave, or (ii) newly formed dust in the ejecta 
heated by a reverse shock traveling backwards (in mass) into the supernova envelope. Our 
observations show that the bulk of the received IR flux doesn't originate in the metal rich 
parts of the ejecta (eg. the weak detection at the center of the ring). Moreover, these parts 
are expanding much more slowly than the outer parts of the H-rich envelope and they have 
not yet reached the ring. As it is most unlikely that dust formed in the tenuous outer H-rich 
layers of the SN interacting with the ring, we believe that the IR emission is produced by 
preexisting dust in the GSM heated by the outer blast wave (scenario 3-i). 

This interaction converts part of the kinetic energy (KE) into radiative output. That 
leads to a caveat in our reasoning: we assume that all the emission from the ER is in the form 
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of IR emission, and we do not consider the possibility that the dust may just be re-radiating 
a fraction of the ambient radiation if the heating mechanism were in part due to radiation 
and not only to collisions. However, it is most likely that the dust may be shock-heated, 
in which case the IR emission cools the hot plasma that gives rise to the X-ray emission as 
well. Indeed, Dwck ct al. (1987) show that IR emission from collisionally heated dust is the 
dominant cooling mechanism of the shocked gas in SNR. Dwek (1987) computes the ratio 
between the cooling function of the gas via gas-grain collisions and that of the gas via atomic 
processes, as a function of the gas temperature only. Dwek et al. (1987) note that, if the dust 
that gives rise to the IR emission occupies the same volume as the X-ray emitting plasma 
(as appears to be the case in view of our Figure 3), this coohng ratio translates into the 
infrared-to- X-ray flux ratio (IRX). These authors show that the observed values of this ratio 
is significantly larger than unity for the 9 remnants that they considered, and concluded that 
IR emission, mostly attributed to gas-grain collisions, is the dominant cooling mechanism in 
these remnants over large periods of their evolutionary lifetime. 

Park et al. (2003, 2004b) report that as the blast wave approaches the dense circumstellar 
material, the contribution from the decelerated slow shock to the observed X-ray emission is 
becoming significant and increasing more rapidly than ever as of 2002 Dec. 31. They fit their 
data with a two-temperature model for that epoch, and find a temperature of kTe = 0.22 
keV and a luminosity L ~ 1.6 x 10^^ erg for the decelerated slow shock component, and 
kTe — 2.44 keV and L ~ 3.7 x 10^^ erg s~^ for the blast wave shock. Combining our data 
with these results leads to values of IRX ~ 6 and IRX ~ 3 respectively. These values 
fall significantly below the theoretically expected ratios, at roughly the same position as the 
Kepler remnant in Figure 1 of Dwek et al. (1987). These authors argue that in Kepler the 
low IRX ratio may mostly reflect the absence of dust in the material heated by the reverse 
shock. Our observations then suggest that little dust was formed in the cjecta during the 
presupernova mass loss phase of the progenitor. We note, however, that the IRX calculations 
were performed for a local interstellar dust abundance and a standard dust-to-mass ratio, 
which might be unapphcable for SN 1987A. 

The A^-band detections of SN 1987A at days 4100 and 4300 suggested that dust was still 
present in the ejecta and the dominant component of the ejecta's bolometric flux. However, 
the ISO observations were shghtly non-stellar strongly suggesting that some of the dust 
emission was coming from the inner ring, presumably heated by the shock front as the high- 
velocity material hits the ER material (Fischera et al. 2002a). EUV radiation is produced 
in the shock, and ionizes the gas upstream of the shock front revealing the structure and 
properties of the outer nebula. Fischera et al. (2002b) estimate that the CS dust is most 
likely silicate-iron or a silicate-graphite mixture, or pure graphite. These types of grains 
have spectral signatures in the A^-band that were not detected in the dust that condensed 
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in the ejecta (Bouchet 1990), most probably in the form of X-type SiC (Amari et al. 1992; 
Lucy et al. 1991). 

After day 530 the dust emission became the dominant cooling mechanism of the ejecta 
of SN 1987A, radiating away the energy from the radioactive nuclides synthesized in the 
explosion. What happened then to this early dust emission? In Fig. 1 a small flux enhance- 
ment with a flux density of 0.32 ± 0.1 mJy (3cr) can be seen in the center of the ring. It 
is likely that this weak feature is the remains of the dust emission from the condensates in 
the ejecta. Although the filter contains several IR lines such as Arlll (A8.99/xm), SIV 
(A10.51/um), and Nell (A12.81/im), these lines require too high ionization state to contribute 
significantly to this emission. Moreover, we note from Fig. 2 that the ejecta has not been 
detected neither in the He I line, nor in the Ha or [N II] lines. As the ejecta has not been 
detected in the radio regime either, the possibility of free-free or synchrotron emission is also 
ruled out. Therefore it is likely that the weak feature detected at 10//m is the remains of 
the dust emission from the condensates in the ejecta. Probably the ER and the ejecta have 
different dust temperatures as was observed in Cas A (Dunne et al. 2003). We do not detect 
the ejecta in our Qa-band image, which sets a flux upper limit of 41 mJy at the 3-0" detection 
level at this wavelength, and then a lower limit for the temperature of the ejecta T > 90 K. 
Fransson & Kozma (1993) showed that time-dependent effects due to long recombination 
and cooling times lead to a frozen-in structure of the ejecta of SN 1987A. However, Frans- 
son & Kozma (2002) claim that freeze-out is important for Hydrogen at times shorter than 
~3000 days, and that at later stages positron input from ^^Ti and reverse shock heating are 
the major sources of energy input to the envelope. In view of their Figure 2, the bolometric 
luminosity of the ejecta should be of the order of 1. — 2. x 10^^ erg s~^. Assuming that there 
is no other energy input than the radioactive decays, this luminosity is an upperlimit of the 
IR luminosity (e.g. all the energy might not be released in the IR through dust emission). 
To derive this luminosity given our 10 fim. flux measurement, the dust temperature must be 
T < 100 K. Thus, our observations give us upper and lower limits, implying a temperature 
90 < T < 100 K for the dust in the ejecta. The dust mass derived from equation (1) is then 
M — 10~^ — 2. X 10~^ considering the same approximate range for K(lOyum) of 100 — 2000 
cm^ g~^ as for the ER. 

We stress that the value given here is an upper hmit of the IR luminosity, as Dwek et 
al. (1992) show that at day 1153 hue emission, mostly due to the [Fell] 26^m transition, and 
other ground-state fine structure lines at 35, 51 and 87 yum, may contribute significantly to the 
total luminosity. These authors find for that epoch a dust temperature 120 < T < 160K. 

We note that the same caveat given for the ER applies to our statement concerning the 
energetics of the ejecta: we do not consider the mechanical energy, which is in the process of 
being transferred to the ambient medium through the forward shock or to the ejecta itself 
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through the reverse shock, that can power the observed IR emission. This issue will be 
tackled in a forthcoming paper. 

Fig. 4 shows that the flux at 10 fj,m dechnes exponentially from day 2200 through day 
4200 (the ISOCAM and OSCIR observations). This figure shows two data points for day 
6067 (ejecta and ER). The possible ring contribution to day 4200 is unknown. However, 
it seems that the ejecta 10 ^m flux from day 6067 lies on the extension of the exponential 
decline curve that roughly goes through the 10/xm flux at day 4200. If one assumes then 
that the lOyum emission from the ejecta indeed follows an exponential decline law, that leaves 
no room for any ring contribution to day 4200. We conclude therefore that the 10 fim. flux 
decline of the ejecta from day 2200 to 6000 is roughly 0.32 mag y^^, and that the ring 
emission started most likely around day 4000 or later, in good agreement with Fig. 5 of Park 
et al. (2002) for ATCA and Chandra/ROSAT data. 

Gerardy et al. (2000) report on the detection of dust emission in the near-IR spectra 
of the Type Iln supernovae SN 1998S (a supernova similar to SN 1987A) and SN 1997ab, 
and summarize dust emission in other supernovae: the IR excess observed in SN 1979C and 
SN 1985L is interpreted as dust forming in the ejecta, while it is attributed to dust lying 
in a preexisting GSM in SN 1982E, SN 1982L, SN 1982R, SN 1993J, SN 1994Y, and SN 
1999el (Di Garlo et al. 2002). Although it has been claimed (Dwek 1983; Dwek et al. 1983) 
that the ncar-lR emission from SN 1980K can be explained as an IR echo, we note that this 
supernova displayed, together with an IR excess, large blueshifts of the [O I] and [O 111] lines 
of ~3300 km s~^ (Fcscn & Backer 1990) compared with 1870 km s~^ in SN 1987A, which 
seems to imply that dust has formed also in the ejecta of SN 1980K. Also, Elmhamdi et 
al. (2003) presented clear evidence for dust formation in the Type IIP SN 1999em 500 days 
after explosion. Lagage et al. (1996); Douvion et al. (1999); Arendt et al. (1999) detected 
a continuum consistent with silicate dust of very small pyroxene grains (MgSiOa) in the 
youngest known Galactic SNR Gas A. 

All of these above observations of dust in supernovae may lead us to the conclusion that 
supernovae could be a major source for interstellar dust. Various indirect arguments (Amari 
et al. 1992) and theoretical calculations (Dwek & Scalo 1980; Glayton 1982; Dwek 1998) 
also suggest that hypothesis. Nevertheless stellar winds from various types of stars are still 
considered a major source of dust production. For instance, Douvion et al. (1999) argue that 
supergiants and AGB stars form the bulk of the dust in the Galaxy, not supernovae. Dunne 
et al. (2003) argue that SNe are at least as important as stellar winds in producing dust in 
our Galaxy and in galaxies at high redshifts. Morgan et al. (2003) also conclude that SNe, 
or their progenitors, may be important dust formation sites. 
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4. Considerations for the Future 

Several theoretical models predicted the presence of dust in the CSE of SN 1987A which 
was produced by a wind in the supergiant phase. The CS dust is likely to be in the form 
of pyroxene or graphite grains, while the dust condensed in the ejecta is most probably of 
some type of silicon carbide. Forthcoming imaging and spectroscopic observations focusing 
on line versus continuum emission should determine the type of emission present, and are 
required to clarify the issue of the energetics of the supernova. 

Finally, in order to assess the role of SNe in the production of dust in the Universe, it 
is clearly important to measure the presence of dust that survives into the formation of the 
remnants, and for this, mid-IR and sub-mm observations are critical. 
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Fig. 1.— SN 1987A seen with T-ReCS at day 6067 in the N band filter (A = 10.36^m). The 
image is smoothed 2 pixels (0.18 arcsec). Note, in particular, the central point source that 
corresponds to the ejecta of the supernova 
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Fig. 2. — T-ReCS lO/iin images of SN 1987A smoothed to 0.18 arcsec (2 pixels) resolutions 
are overlayed with contours of the images obtained in the Hel line (1.083/im) with OSIRIS 
at the CTIO Blanco telescope at day 5749 (a), and through the F656N filter (Ha and [N II] 
A6583 A) with HST at day 5555 (b). The HST image is smoothed 1.5 pixels (0.15 arcsec). 
The central source is not seen in any of the two overlays while it is detected at 10/xm. 
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Fig. 3. — SN 1987A: the Chandra (a) 0.3-8.0 keV broadband image obtained at day 5791, 
and the ATNF 8 GHz restored image (b) obtained at day 6002. The overlays correspond to 
the T-ReCS image. 
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Fig. 4. — Mid-IR light curves of SN 1987A. The point called 'Ejecta' derives from the weak 
point source near the center of the image 



